The spectral dimension of random trees 
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We present a simple yet rigorous approach to the determination of the spectral dimension of 
random trees, based on the study of the massless limit of the Gaussian model on such trees. As 
a byproduct, we obtain evidence in favor of a new scaling hypothesis for the Gaussian model on 
generic bounded graphs and in favor of a previously conjectured exact relation between spectral and 
connectivity dimensions on more general tree-like structures. 
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INTRODUCTION AND SUMMARY 



The spectral dimension d was first introduced by Alexander and Orbach [Qj to characterize the low frequency 
vibration spectrum of fractals. Nowadays, it is considered by most the right generalization of the euclidean dimension 
of regular lattices to irregular structures in general, whether fractal or not, such as polymers, glasses, percolation 
clusters, dendritic growths and so on. This is verified at least in many experimental situations, as well as in some 
theoretical model examples, whenever the physical dimensionality explicitly enter observables such as the spectrum of 
^ ■ density fluctuations, the long time properties of random walks, the critical exponents of statistical models and many 
CO | others. 

In a theoretical setting, we may consider a connected graph as model of generic irregular structure, with physical 
microscopic degrees of freedom attached to the nodes and interaction among them associated to the links. This 
scheme implies that the coordination of nodes is bounded (the number of nearest neighbors of an atom, molecule or 
' basic building block has a geometrical upper bound) and, if the graph in question is infinite, that the surface of any 
Van Hove ball (defined in terms of the chemical distance alone) must be negligible with respect to the volume when 
■ its radius goes to infinity. In other, more physically immediate words, we say that the structure is embeddable as a 
whole in a finite-dimensional euclidean space and a standard thermodynamic limit of infinite systems exists. We are 
thus led to consider bounded graphs || and to study their spectral dimension. 

To be precise, we refer here to the average spectral dimension, that is a global property of the graph, related for 
instance to the infrared singularity of the Gaussian process or, equivalently, to the density of small eigenvalues of the 
Laplacian or the graph average of the long-time tail of the random walk return probability On macroscopically 
£j , inhomogeneous structures this average dimension may differ from the local one, which enters the long-time tail of 
the random walk return probability on a given node []3|. However, it is the average spectral dimension that plays the 
same role of the lattice euclidean dimension in many contexts, as for instance in providing a consistent criterion on 
whether a continuous symmetry breaks down at low temperatures . 

In this work we address the problem of determining the average spectral dimension of bounded trees and in 
particular of statistically homogeneous random trees |2j (this homogeneity actually implies that local and average 
spectral dimensions coincide), which are defined by the following branching algorithm: given a set of coordination 
fractions {/ z , 1 < z < z} that are properly normalized and have average coordination 2 

z z 

£/, = l , E z ^ = 2 (!) 

z=l z=l 

an infinite random tree is built by first selecting with probability f z the coordination of the origin and then, proceeding 
shell by shell away from the origin, by choosing the coordination of every new node with probability f z except for 
the first node of each shell, when the modified probability f z = (z — l)f z is used to implement the non-extinction 
preconditioning (this is algorithm B of pi). These random trees represent a generalization of the incipient infinite 
percolation clusters on Bethe lattices || , since the branching probabilities are freely assigned rather than being fixed 
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by a single percolation probability. They also generalize to arbitrary coordination fractions f z the infinitely large 
specimens of the branched polymers studied for instance in . 

Since the introduction of spectral dimension, the value of d = 4/3 has been conjectured for the incipient infinite 
percolation cluster in any dimension. This conjecture was proved false at low dimensions, while has been considered 
valid in the limit of infinite dimension (percolation on Bethe lattices). Leyvraz and Stanley had in fact produced 
an heuristic scaling argument in favor of d = 4/3 already in 1983. Several years later, a more rigorous approach 
was pursued in where the branched polymer phase of two-dimensional quantum gravity was proved to have 
indeed d — 4/3. Branched polymers correspond to a grandcanonical ensemble of trees with all possible sizes, which 
is generated by a branching process like ours above but without the non-extinction preconditioning. In this case 
the definition itself of spectral dimension is subtle since such concept does not apply to finite graphs. This is the 
major difficulty in the approach of ||. Another mathematically well founded approach to the problem is given in 
where the exponent 9 of the anomalous diffusion on our class of random trees is determined as 9 = 1 (the author 
reports only the shorter and simpler of two different derivations, which turns out to be more than twenty pages long 
and states that the other, more complete treatment is "monstrously long"). It is commonly held, although never 
rigorously established in general, that 9 is connected to the spectral dimension and the (average) connectivity (or 
intrinsic Hausdorff) dimension d c as: 

1 1 (2) 



2 + 6> 2 d c 

Therefore 9 — 1 would imply d = 4/3, since d c = 2 for random trees Q. 

In this paper we provide a new, independent and fairly ri gor ous derivation of d — 4/3 based on the infrared 
properties of the Gaussian model on our random trees (section |Vj) . The connection of the Gaussian model with the 
random walk on any graph is a well known fact that we briefly review for completeness in section |b|. As a byproduct, 
we obtain also an argument in favor of the conjecture (||), for any infinite graph for which the thermodynamic limit 
of the Gaussian model exists: it is natural in fact to expect the scaling form 

li ~ /ior d "<=/(/i r 7 ) (3) 

for the effective squared mass of the Gaussian model in terms of a very small "bare" squared mass /xq and the very 
large radius of a graph (see section [b] for the proper definitions); moreover we conjecture also that the finite-size 
exponent 7 coincides with 2 + 9, which is called sometimes random walk dimension and characterizes the scaling 
t ~ r 2+e between the time needed by a random walker to visit with non-negligible frequency nodes at distance r from 
its starting point; then the existence of the thermodynamic limit on the effective mass implies that f(x) ~ x^ dc ^ 
and so 

M ~ Mo 

yielding exactly equation (0), since in general fi ~ //J d ^ 2 as /io — > (see section ||). 

In our treatment of random trees we are indeed able to explicitly derive the scaling law (|^) with d c = 2 and 
7 = d c + 1 = 3, demonstrating d = 4/3. It is interesting to notice that if the scaling law (m is assumed a priori, with 
d c = 2 as well known for random trees, then the precise numerical value of 7 (and therefore of d) follows from a very 
simple "geometrical" fact: as shown in section |v|, the first correction to the finite tree behaviour fi oc /ior^ c comes 
from the sum of squared volumes of all subtrees, which grows like r 2dc+1 if the volume of tree grows like r d *. Then 

d c + 1 

which is also a known conjecture, even with d c =/= 2, for the incipient percolation cluster and other tree-like disordered 
structures ||. Notice that d is monotonically increasing from 1 to 4/3 as d c grows from 1 to 2. 

Let us observe that d c = 2 and d = 4/3 are the greatest value among all known average connectivity and spectral 
dimensions for bounded trees. This fact lead us to formulate the following conjecture: for all bounded trees d c < 2 and 
d < |j with saturation only in the case of random trees. Evidently this must be true only for the average dimensions; 
there indeed exist examples of trees (such as NTp [|Io[ ) whose local connectivity and spectral dimension, d c and d, 
exceed 2 and 4/3, respectively; these present however macroscopic inhomogeneities so that the local dimensions are 
different from the average ones which do satisfy the bounds at 2 and 4/3. Let us observe also that equation ([|) does 
not apply in the same examples if the average dimensions are replaced by the local ones. 

At the moment these upper bounds are only conjectured, but it has been rigorously proved that d < 2, since 
bounded trees are always recurrent on the average (TTJ . 
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II. GAUSSIAN MODEL, RANDOM WALKS AND SPECTRAL DIMENSION 



The Gaussian model on a generic connected infinite graph Q is dehncd |12 by assigning a real- valued random 
variable <f> x to each node x € Q with the following probability distribution 



dv r [<t>] = J- CX P ~\ E <f>x(L xy +Mo^)^] Yi d< ^ x ( 5 ) 

x,y£B(o,r) x^B(x,r) 

for the collection <j> = {<p x ; x G B(o, r)} relative to any Van Hove ball 

B(o, r) = {y € Q : d(o, y) < r} , o E Q , < r < oo 

where d{o, y) is the chemical distance. In equation (|^) Z r is the proper normalization factor (the partition function) , 
fio > is a free parameter (the squared mass in field-theoretic sense) and 

-Lxy — z x S xy A X y (6) 

is (minus) the Laplacian matrix on Q, z x being the coordination or degree of x and A xy the adjacency matrix). 

The thermodynamic limit is achieved by letting r — ► oo and defines a Gaussian measure for the whole field <f> — 
{4> x ; x £ Q} which does not depend on the center o of the ball |D| if Q is bounded. The covariance of this Gaussian 
process reads 

{<f) X <f>y) = C X y(Ho) ={L+ fl ) X y (7) 

and setting 

fl - W) _1 A 

Cxy= [ - >XV , W X y = ^^ (8) 

z x + Mo Zy + no 

one obtains the standard connection with the random walks over Q fl^| : 

oo 

a-wOxJ = E( w % = E w m 

t— 7:;e< — y 

where the last sum runs over all paths from y to x, each weighted by the product along the path of the one-step jump 
probabilities in W: 

1 = (x,yt-i,...,y 2 ,yi,y) W[i\ = W xyt _^W yt _ xyt _„ . . . , W y2yi W yiy (10) 

Notice that, as long as (j,q > 0, we have YlxO^^^y < 1 f° r an y name ly the walker has a nonzero death probability. 
This implies that C xy is a smooth functions of fiQ for /io > e > 0. In the limit fiQ —> the walker never dies and 
the standard random walk is recovered; then the sum over paths in equation (^) is dominated by the infinitely long 
paths which sample the large scale structure of the entire graph ("large scale" refers here to the metric induced by 
the chemical distance alone). This typically reflects itself into a singularity of C xy at /iq — 0. Consider for instance a 
diagonal element C xx of the covariance and suppose that its singularity is power-like 

SmgC xx (fi a ) oc ^o /2_1 (11) 

Then by equation (||) one sees that the same power-like behaviour appears, as I — > 1, in the discrete Laplace transform 
Pxx{l) of the random walks return probability P xx (t) at x. In turn this implies the long time behaviour 

P xx (t) oc£- d/2 , t ► oo (12) 

identifying d with the local spectral dimension, which does not depend on the specific node x 113]. To be precise the 
denomination "spectral dimension" refers to the behaviour of the spectral density p(l) of low-lying eigenvalues of the 
Laplacian L, namely p(l) ~ I d / 2 ^ 1 . In fact, it can be shown (j| that d is connected to the long time behaviour of 
the graph average of P xx (t), or equivalently to the infrared power-like singularity of the graph average C(/iq) of C xx , 
that is 

Sing C(^= Sing lim 1 V C**(mo) cx (mo)^ 2 " 1 (13) 
r-+oo \B(o. r) 

1 V y x£B(o,r) 
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where |S(o,r)| is the volume of B(o,r). d is therefore called sometimes "average" spectral dimension. On regular 
lattices it coincides with the local spectral dimension and the usual Euclidean dimension. 
By definition, the diagonal elements C xx of the covariance may be written as 

C xx = f d<p<p 2 Z x (ip) (14) 

where 

Z x (ip) = lim / dv r [4>] 8(4> x - ip) 



Thanks to the fundamental self-reproducing property of Gaussian integrations, Z x (ip) is necessarily a normalized 
Gaussian in (p, that is 



■<*>-(£ 



1/2 , , 



where, for consistency with equation (14) 



is the "effective mass" at x. 



» M = ^M (15) 



III. THE GAUSSIAN MODEL ON RANDOM TREES 



Let us now consider the Gaussian model on an infinite bounded tree T. Any node x € T may be regarded as the 
root of an infinite branching process which produces T as family tree rooted at x. This simply means that x has as 
many descendants, or children, as its coordination z x , while any other node x' has z x i — 1 children one step further 
away from the unique ancestor x and a father one step closer to x. In particular, the z x children y of x are themselves 
roots of disjoint subtrees, or branches T', which may or may not be infinite. Now, by the simple rules of Gaussian 
integration, one easily verifies the following relation between the effective mass at x, relative to the entire T and the 
effective masses at children y of x, each one relative to the corresponding branch Ty : 

a children +M ^ [ y J 
of x 

A similar relation holds for the effective masses ^ y [Ty] in terms of the z y — 1 masses of y's children, and so on for the 
rest of the tree. 

Next, suppose that the branching process is the Galton-Watson random one, with the non-extinction precondi- 
tioning, described in the Introduction (and in more detail in ||), so that T is an infinite random tree. Because of the 
preconditioning, T may be regarded as the union of infinitely many subtrees generated by independent and identically 
distributed (i.i.d.) random branching processes, without the non-extinction preconditioning, whose roots are attached 
to the nodes of a half-infinite chain (the spine) that starts from x. This means that only one of the branches at a; is a 
priori guaranteed to extend to infinity, while the other z x — 1 are the family trees of i.i.d. random branching processes 
which might stop after finitely many generations. 

Let us consider first the probability distributions relative the branching process itself and in particular the probabil- 
ity Pr(p-) that the effective mass at the root x be exactly fi when the branching algorithm is iterated for r generations, 
that is 

P r {n) = {6{nx - n)) r 

where the average is over all realizations of the tree with r generations, namely on all possible histories of the branching 
algorithm with r shells completed. Taking equation ( |l6| ) and the spinal decomposition into account, we may write 

/z z 
dmdv* ...dM(/*-w>-£ Y^-) p 'r^) II p ° (Mi) 

3=1 ''' J=2 
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where P°(fi) is the same probability without the non-extinction preconditioning and fulfills the functional recursion 

K+M = £/» / ■ ■ ■ d/*»-l " MO - Yl ) II P °(W) ( 18 ) 

2 J 3=1 3=1 

while P*(/z) is the probability when also the coordination of the root is extracted with the modified probability f z 
and satisfies therefore 

« z — 1 z—1 

P'r+M = £/« / • ■ • d/i*-l *(a* " A*0 " £ TT^) P '(Ml) II Pr ^) (19) 

2 ■/ 3=1 + M3 J=2 

As r — » oo we expect all these probability distributions to converge to their limiting forms Poo (fi) , P£, (/x) and P^ (/i) , 
turning eqs. jl7|), ( |l8| ) and ([l^) into highly nontrivial coupled integro-functional equations for P 0o (fi) ) P^(/i) and 
P^(/i) with a parametric dependence on /iq. In particular we are interested in their infrared behaviour as fio — > 0. 



IV. STATISTICAL HOMOGENEITY AND AUTO-AVERAGING PROPERTY 

Before addressing the just mentioned proble m, we need to establish the precise relation with the original problem, 
that of determining the fio — ► behaviour of C(/io) on a given "generic" infinite random tree T produced by our 
algorithm. By "generic" in this context we mean an infinite T belonging to the subset of unit measure of trees that 
fulfill the auto-averaging property for local observables. To be explicit, consider the frequency F(r, r) of occurrence 
in T of a given rooted tree t with r generations. In terms of the Van Hove balls B(x,r), regarded as subtrees of T, 
we may write F(r, r) as the graph average 

F(r, t) = lim F R (r, r) , F R (r, r) = / £ S[B(x, r) = r] 

where o is an arbitrary node of T and S[t' = r] is one if its argument is a true statement and zero otherwise. F(r, r) 
does not depend on the choice of o (T is a bounded graph) and it is convenient to identify it o with the root, or origin, 
of the infinite branching algorithm. 

Let us then recall that any T produced by the algorithm may be decomposed into the spine S (the origin o plus 
all the nodes on which the coordination was extracted with f z ) and the transverse family trees of the identical, 
independent and non preconditioned branching processes starting from the nodes of the spine (see Q || for more 
details). Evidently, the interesection of the transverse family trees rooted at s 6 S with B(o,R) contains at most 
R — d(s,o) generations, so that they contribute differently, even on average, to Fr(t,t). Therefore, to simplify our 
derivation, we consider a modified thermodynamic limit, in which the infinite tree is not recovered through a Van 
Hove ball with diverging radius, but as R — > oo limit of the set 

= [J W S ,R 

S&Sr 

where Sr is the restriction of the spine to the first R generations of T and lj Si r are the union of all transverse family 
trees rooted at s and restricted to the first R generations. Thus we redefine Ffj(r, r) as 



Pfl(r, t) = "52 w s ,r Q s ,fl(r, t) , w S! r= *' R , Q s .r{t,t) = V 5[B(x,r)=T] 

SGDH ill x £uj 3t R 



Notice that, by construction, the random weights w s> r are all identically distributed and are all independent except 
for the normalization constraint X) s es R w s.R = !• The random frequencies Q s ^(r, r) are all identically distributed, 
but are not independent because a ball B{x, r) will in general intersect several u) Si r if d(x, s) < r. Therefore we 
split u> s r into ctJ^R, formed by nodes x such that d(x,s) > r, and its complement uf, which does not depend on 
R for large R. This implies an analogous splitting for the random weights and the random frequencies. Now the 
random frequencies R (r, t) are all independent since the two events B{x, r) = r and B(x' ', r) — r are independent 
if x G and x' S u>^, R with s ^ s' . Notice also that Q s ^R(r,r) and w s > t R are independent when s ^ s'. 
Thus we have 

F R (r 7 r)^F^(r,r) + F>(r,r) 
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with Fft (r, t) getting contributions solely from nodes at distance not larger than r from the spine. Clearly, in the 
limit R — > oo, these nodes form a subset of zero measure of 7", since lj s ojf{r, t) is essentially one dimensional while 
a generic T has d c = 2. Hence, as R — > oo, 

F(r,r) = lim F>(t,t) 

R — >oo 

By analogous dimensionality arguments we get that almost always liiriR_ >00 w a r = 0, although the sum of all these 
random weights stays properly normalized. 

We may now invoke directly the law of large numbers for Fo{r, r), which is the weighted average of i.i.d. random 
variables. Hence, for any finite r, F(r,r) is almost always non-fluctuating on infinite trees and 

F(r,T) = (Q>(r,T))° 

where Q > (r, r) is anyone of the liniR-^oo R (r, r) and the average (■) is taken with respect to the algorithm without 
the preconditioning on non-extinction. But any node x on such trees laying at a distance larger then r from the origin 
enjoys statistical homogeneity for the event B(x,r) = r, since the branching is i.i.d. on every node (but the origin, at 
most) and B(x, r) = r cannot contain the origin. Moreover, the unique path connecting x to o plays the role of the 
spine for a branching process rooted at x and preconditioned on non-extinction for at least r generations, since the 
trees contributing to (Q > (r, t)) have at least r generations. We conclude therefore that almost always an infinite 
random tree T has the auto-averaging property 

F(r,r) = (d[B(o,r)=r]) (20) 

with respect to the probability itself that the first r shells of T build up to form the given finite tree r. Then 
this property holds for any local observable and in particular for the effective mass fi x [B(x,r)] of a Gaussian models 
restricted to B(x, r). Thus the probability distribution P r {^) coincides with the frequency of the event [i x [B(x, r)] = /(. 
as x varies throughout a "generic" infinite random tree T. The same would evidently apply to the probability 
distribution of the variable (2/i) -1 and the frequency of the event C xx [B(x,r)} = (2/i) -1 . 



V. SCALING 

Let us now return to the recursion relations (|l7|), (18) and ( |l9| ) for the probability distributions P r (/x), -P°(m) and 
P*(fj,), respectively. We are interested in the thermodynamic limit r — > oo of infinite trees followed by the massless 
limit /Lto — > 0. This two limits can actually be applied simultaneously, provided we identify the correct scaling variables. 
Suppose in fact that, as r — > oo and fi, — > 0, the scaling law (|3|) applies; since fi is a really a random variable, a 
better formulation is in terms of its probability, that is 

P r (n;no) - Mo ^(MMo,Mor 7 ) (21) 

for suitable scaling exponents (3 and 7 and scaling function F(u, s). (3 is then directly related to the spectral dimension, 
as (3 = —l + d/2, thanks to equation (|ll|), the results of the previous section and the existence of the thermodynamic 
limit; 7 provides instead a measure of finite size effects at very large sizes. In particular, if we let fj,Q ~ > too fast, so 
that /ior 7 — ► even if r — > 00, we must recover the massless behaviour characteristic of finite trees. This amounts to 
linearize the effective mass formula (|l6|), so that 



where v[T] is the volume of T. Thus 



Hx[T\ = hqv[T] 



Pr(MIM)) ^ Mo ^r(/t/W)) = Mo lp r{v) 



where P r (v) is the volume probability at radius r, which for large r has a scaling form in terms of v/r 2 explicitly 
determined in (thus implying d c — 2). Since this scaling form is highly nontrivial, with exponentially small 
behaviour also for v/r 2 — > 0, compatibility with the scaling hypothesis ( ^l|) requires 1 + /3 = 2/7, that is 

1 _ d _ d 



7 



One is therefore lead to the identification (compare with equation (g) in the introduction) 

7 = 2 + 

relating finite-size effects to the exponent of anomalous diffusion. Of course, since the relation (^) is only a well 
founded conjecture, an explicit derivation of 7 = 2 + 8 would provide a proof of the conjecture for our random trees. 
Let us stress the point: 7 measures the scaling behaviour of the Gaussian model (and therefore of the random walk) on 
finite random trees of very large radius r; 2 + 6 characterizes instead the scaling t ~ r 2+e between the time needed by 
a random walker to visit with non-negligible frequency nodes at distance r from its starting point. The two concepts 
are certainly related but quite distinct, so that the identification 7 = 2 + has a nontrivial content. 

Let us come back now to the opposite limiting behaviour of the probability distributions that is the subject of this 
work, namely r — * 00 and /£, /l«o — » with u = (j,/j,q fixed and s = /iof 7 finite and possibly very large. We shall show 
now that the scaling law (|2l| ) is indeed correct with the precise values 

(3 = -1/3, 7 = 3 (22) 

implying d — 4/3. The detailed proof, by induction on all moments of P r (/i), P°(p) and P'(p) is rather simple albeit 
notationally a bit cumbersome and is reported in the Appendix. Here we only provide a sketch of the derivation. 

The first step consists in rewriting the recursion relations (0), ( |l8| ) and ( |l9| ) into a more tractable form. In fact 
they are amenable to a direct treatment via Laplace transform only in the limit fiQ — > at fixed r, while we are 
interested in a regime where the nonlinearity in the effective mass formula (|l6|) plays a crucial role. Therefore we 
consider the expansion in all powers of fio 

00 

^[T]=^(-ir +1 K(xK ( 23 ) 

71=1 

where the coefficients V n (x) are integer-valued random variables satisfying composition laws obtained by equating 
powers of ^0 on both sides of equation (|l^). These read 



vi( x ) = i+ J2 v *(y) 



y children 
of x 



y children 
of x 



Va(x)= [Vz{y) + W 1 {y)V 2 {y) + V 1 {yf] 



y children 
of x 



v *l*) = E [ y 4^) + 2V ^(y) v 3(y) + y^ivf + ZVx(y) 2 V2{y) + Vi (y) 4 



y children 
of X 



V 5 (x) -- 

which we write in general as 



V n (x) = 5 nA + ]T F n (Vi(y),V 2 (y),...,V n (y)) (24) 



y children 

of x 



where the form of T n can be easily induced from the previous formulas. Notice that V\ can be interpreted as the 
volume of the subtree, Vi as the sum of the squares of all sub- volumes and so on. Now we may write recursion 
relations for the probability of the coefficients V n {x) themselves. For instance, in the branching algorithm without 
non-extinction preconditioning, we have 



2—1 n z—1 



P° +1 (V 1 ,...,V n ) = J2f z ]T ]JP°(V 1 , i ,...,V n , i )l[S(V k -5 ntl -J2MVi,h--->Vk, j )) (25) 



z Vi 1 — Vi,»_i i=l k=l j=l 

V 2 ,i-V 2 ,z-i 

v„.i"-i>„.*_i 



with similar expressions for probabilities P r (Vi, . . . , V n ) and P*(Vx, . . . , V n ) corresponding to those for the effective 
mass Pr(lJ-) and P*(/x). For brevity we do not write them here; they can be found in the Appendix. 
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The advantage of this new formulation stands in the fact that T n {{V\^, . . . , V n ,i) is by construction linear in V n ^. 
This allows to set up a systematic control of the large r behaviour of the recursion relations for the multiple moments 
of P r (Vi, . . . , V n ), P°(V%, . . . ,V n ) and P*(V%, . . . , V n ). Consider for example the lowest moments: from equation (pH) 
and its analogs in the Appendix we get 

<^>; +1 =e^( i +(- i )(^>:)= i +(^>; 

z 
z 

= e /. f 1 + w; + (* - D<%>;) = 1 + w; + wi 



(26) 



(27) 



where g(Z) is the so-called probability generating function 

z 

which thanks to eqs. (Q) enjoys the properties 

5(1) = 1, g'(l) = l 
Hence we obtain the asymptotic behaviour for large r 

(^»r, ( Vl ) r ~ r 2 (28) 

which implies d c — 2 since Vi is just the volume. Next consider the averages of V2 and V] 2 ; the recursion rules read 

z 
z 

= (V 2 )' r + (V?) m r + g"(l)(V 2 )° r + g"{l){V? )° r 

<^ +1 =<^>;+<^>;+<^>;+<^ 2 >; 

and 

(^ 2 >: +1 =i+^i 2 >;+5"(i)(^i>;) 2 +2<^>; 

<v?>; +1 = 1 + (v^: + + /a) [<^ 2 >: + + w;] + axw;) 8 
= 1 + (v?)' r + 2( Vl ): + ( Vl y r + + 2{ Vl y r + g'^my 

We can read right away the r — > 00 asymptotic behaviour 



(29) 
(30) 



Evidently the moments of P'{V\, ... ,V n ) have all the same asymptotic behaviour of those of P r {V\, ... ,V n ) and need 
not be written out explicitly. 

Let us make use of the results (^8|) and ((2^) for the mean value of the effective mass at the root o: 



<M> r = E(- 1 )" +1 ( y «)^o - Vor 2 (c + c W0 r 3 + . . .) 

n=l 

where cq and c\ are constants. If we assume the scaling law (|l|) for the probability P r {^', Mo), then clearly 



(31) 
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with Fi(s) finite in the limit s — > oo; consistency with ( |3l| ) now requires that F\(s) ~ s 2 / 7 (l + 0(s 3 )) as s — > to 
reproduce the correct r behaviour on large but finite trees; but matching also the powers of /io implies 

1 + = 2/ 7 , 2 + p = 5/ 7 

yielding exactly the desired result (|2^). This argument may be turned around to avoid referring to the scaling 
hypothesis (|2l]) for the probability and to better clarify the role of the thermodynamic limit. In fact, the first two 
terms of the /io— expansion in (|3l| ) suggest that, as r — > oo and /io — > at fixed /io^ 3 > 

(M> r ^Mor 2 Gi(Mor 3 ) (32) 

for a suitable scaling function Gi(s). Now it is the existence of the thermodynamic limit [|12| requires that Gi(s) ~ 
s -2/3 f or i ar g e Sj yielding again (fi) ~ /i^ 3 , that it /3 = —1/3. Of course this formulation is only tentative, since we 
still need to verify that, for any other n > 2, the n— th term in the /io— expansion of (/i) contains r 3rl+2 as highest 
power of r. Moreover, even if we can prove the scaling J32] ) of (/i) , this does not imply the scaling ([H]) for the whole 
probability Pril^l A*o)- I n the Appendix we demonstrate by induction that, as r — > oo 

k 

<n^)r~ r8Ein '" fc [ i +°( r ~ 1 )] ( 33 ) 



(M) r - cofJ, r 2 + ci/^r 5 + c 2 ^lr s + ... + Ci/j^'r 



Hence, by retaining only the dominant r— behaviour, we have 

- 2 + Cl rf~ 5 — 3 - 8 ^ 1 

~ /i r 2 Gi(/i r 3 ) 

and, by the same token 



oo k— 1 

= E(-!)^o E <^^-i) r - Mc> 4 G 2 (/ir 3 ) 

fc=2 j=l 



Similarly one obtain for any n 



(n n ) r ~ ^r 2 »G„( Mo r 3 ) 
so that, requiring the existence of thermodynamic limit, one obtains as /io — > + at r 

( M «) = lim (/i") r = fe„/i l/3 

where 



We may also reformulate this result as 



lim s 2 "/ 3 G n (s) 



</i«) r ^/i "/ 3 ^(/i r 3 ) 



where F n (s) = s 2n ^ 3 G n (s) has a finite nonzero limit as s — > oo. This general form of the moments of P r {n', Mo) implies 
the scaling form ( pi] ) with (3 — — 1/3, 7 = 3 and 

F„(s) = J duu n F(u,s) 

for the probability -P r (/i; /io) itself. 

Let us close this section with one more observation concerning the universality of the scaling law (||) or (|2l]) with 
respect to the coordination distribution defined by the f z numbers or, equivalently, by the probability generating 
function g(l). A more careful analysis of the dominant large r terms in the recursion rules for V n expectations shows 
that (see Appendix for details) 

k 
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FIG. 1: Lines of constant integrated probability from numerical data. 



where a = g"(l) measures the fluctuations of the coordinations since 

a = ]T(z-2) 2 / 2 

z 

We may therefore refine the scaling law (|l]) to 

PrO;Mo) ^ (aMo) _1/3 i ? ((aMo)" 1/3 M I a/i r 3 ) 
where F (u, s) is a universal scaling function. 



(34) 



VI. NUMERICAL CHECKS 



As an extra check on our derivation and to determine the profile of the universal function F(u, s), we numerically 
evaluated P r (fJ,', no) by directly calculating the effective mass // for each member of a large set of trees produced by 
our branching algorithm. For each one of the selected values of no and r we obtained several thousand values of /i 
that may be arranged into an histogram providing a discretized estimate of P r (n] no)d-n- This procedure does not 
require much computer time or memory thanks to the recursive structure of the mass composition formula (|l6|). Our 
results for the integrated probability at the five values 0.05, 0.2, 0.5, 0.8 and 0.95 are plotted in figure pi using the 
scaling variables u = /i /i and s — no 1 " 3 (with our specific choice of f z we have a — 1). The data show a very good 
scaling profile with a very well defined crossover, as s gets larger, from the behaviour typical of small trees, n oc no, 
to that characteristic of large trees, that is /i oc /i^ 3 . The vertical section of data at s ~ 10 4 provides a good estimate 
of the probability distribution F(u, oo)du and is reported in figure ^| 
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FIG. 2: Numerical distribution for F(u,s) for 10000 random trees. 



VII. CONCLUSIONS AND OUTLOOK 



Let us remark once more on the two distinct levels of our approach: our proof by induction of equation ( p3| ) provides a 
rigorous basis for the scaling form (21 ) of P r (/x; hq) with (3 = — 1/3 and 7 = 3. Thanks to the auto-averaging property, 
this implies in turn that indeed 

Sing C(/i ) oc Mo 1 i d = 4/3 

Moreover, due to the statistical homogeneity of our random trees, we see that the local and average spectral dimension 
coincide. On the other hand, if we assume scaling, as natural and common in this contexts, then the determination 
d — 4/3 follows solely and very simply from the large r asymptotic behaviour of {V2) , that is the average of the sum 

of squared volumes of all subtrees of the tree. Since the volume Vi scales as r d °, with d c = 2, then V2 must scale, 
almost by definition, as r 2dc+1 ; then 



{fi) r ~ /lo^ (c + c lf i r d ^ +1 + ...) 



which is consistent with scaling only if 



d = 



2d r 



1 



(35) 



Notice that this holds true even if d c < 2, as could happen for different types of bounded trees. We conjectured in 
Pj that d c < 2 for any bounded tree, with saturation only for random trees. Then automatically d < 4/3 for any 
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bounded tree, with saturation only for random trees. The general expression ( |35| ) was conjectured years ago from 
very different considerations || . 

Another relevant issue concern the spectral weight w of infinite random trees. This is the coefficient of the 1 
singularity of C(/xo) as — ► 0, that is 

C{no) ~ w^ 1/3 

It is clear from equation ( |l5|) and the auto-averaging property that the spectral weight is related to the scaling 
function F(u, s) by 

/du 
—F(u : oo) 

On the other hand it is known p"2| that on any given infinite graph with local spectral dimension d < 2, the local 
spectral weight does not really depend on x, that is 

C xx {fi ) ~ w^t 12 ^ 1 ( 36 ) 

for any node x. Moreover, thanks to statistical homogeneity, it is clear that on a given tree w = (2-u) -1 , with u a 
specific instance of a random variable distributed according to F(u, oo). Since F(u, oo) is not 5— like (see previous 
section), we have in general w ^ w, that is local and average spectral weights differ in general, unlike the spectral 
dimension. One might wonder then about the auto-averaging property, since w does not fluctuate over the tree. Of 
course the resolution of this apparent contradiction is in the order of limits: averaging over the infinite graph and 
small /io asymptotics of C xx (hq) are operations that cannot be commuted in general. 

These last remarks suggests a further development in the explicit calculation of the universal limit function F(u, oo). 
Another interesting direction of research is in the construction of modified tree-generating algorithms such that scaling 
is preserved with d c < 2, in order to check the general formula (|35|). 
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APPENDIX A: ASYMPTOTIC BEHAVIOUR OF MOMENTS 

First of all let us write the explicit form of T n 

n n—1 n—X k k 

Fn(v 1 ,v 2 ,...,v n )=v n +Y^ e ■•• e <^-x>«)n y «i 

k—2n\ — l rik — 1 j—1 j—1 

In the case without the non-extinction condition we will prove that 

< n v ^)i - aY " tij_i r3Zj nj ' h ~ i ( A1 ) 



First of all let us prove this equation when all n 3 ■ — 1; the case k = 1 is known from equation ( |28| ) while for arbitrary 
k is proved by induction; in fact if 

{Vl)° r ~ a^- 1 (A2) 
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is assumed for j < k we have 



<^ +1 >: +1 = E/* E 

* Vi,l,Vl t 2 ! ...,Vi,,- l 

= £/. E 

z V r i, 1 ,Vi,2,...,Vi ._i 



2-1 



fe+1 



0-1 



i + e Vi <* n p °(^) 

i=l / i=l 

(£*.<r+*(5> " 



W + 



i = l 



(A3) 



=(r);+«E^>:r H ); 

3=1 

e <^ i >:<^ 2 >:<^ 3 >^+ i -ji-i2-i3)+... 

31 ,32 ,33 = 1 
fe 

+ k(V^)° r + a ^ + lower order terms 

3=1 

where all addends which consist of a product of three or more aver ages and the ones that come from lower powers of 
(53 Vi t i) are negligible in the limit r — > oo thanks to equation (K3), so that 



<^ +i >; 



a (fe+l)-l r 2(fe+l)-l 



extending equation (A2) from fe to k + 1. 

Now, let us prove equation (Al) by induction; in this case the induction is on the quantity TV = 3 5Z^=i n j ~~ ^. Let 
us denote with n the greatest among m, n^, ■ • • , Tik', we have 

k 2—1 z — 1 

< n ^>°+i = e /- e n («w + E^( y M, • • • v nj ,i)) n • • • ^ 



3=1 



2 Vi,i'"Vi lS _i j=l i=l 

v a ,i---y ft ,,_i 

fc fe n j - 1 

(n^);+EE(^vn^);+ 

3=1 3=1 ' = 1 

"E (II Vn i)°r( II Vn iTr + l0Wer OTder tei " mS 

s jes as 

k 
3=1 



where 5 are all the possible ordered subsets of {ni,ri2, . ■ . ,rik}- Now equation (Al) follows but a few remarks are 
needed. First, all terms in the right hand side involve averages whose value of N is smaller than the one in the left 
hand side. Then all the terms that are not explicitly written are proportional to lower power of r: some of them in 
fact consist in the product of three or more averages (each of them carries a "—1" in the exponent), some other are 
expressions with more than k + 1 factors V n , and finally some others come from products involving the S n , t i. 
We turn now to the probability P'(Vi, . . . , V n ) with non-extinction condition for which the recursions reads 



2-1 



2-1 



P; +1 {V X , ...,V n )=Y J h E Wm, • • • . V n ,i) J] Wm, ■ • • , V n ,i) S(v k - S n>1 - £ MVi,j, v K 



Z Vl.l-Vl,,-! 

V 2 ,i-V 2 , z -i 



i=2 



fc=l 



3 = 1 



We will prove that 



3=1 



(A4) 
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Some special cases, including k = n% = 1 have been already shown in section [v|, while the general case is proved by 
induction. In fact the recurrence for ( ^n 3 -) r rea ds (writing only leading order terms, as before) 

( II V^Yr+X = £(* ~ £ II (<W +E^n J (V lli , • • • V^ij) 

3=1 2 Vi,i~Vi,*-i J- 1 i=1 

V2.1---V2.z-1 

Vfi.i-'-V'fi^-i 

2-1 

• • • V n ,l) J] P°(Vl,i, • • • V B ,i) 

i=2 

k k nj-l k 

=<nu+EE <^-« n y v>: + n ^» 

i=i i=i 1=1 j'^j j=i 

aE ( TT TT + l° wer or der terms 

i=i 



E. 3 £ n^— A;— 1 



with the same remarks as in the previous case. 

Finally the recursion for the probability P(Vx, . . . ,V n ) reads 

z n 2 

Pr+x(Vi, ...,V n ) = Y,f* £ P " (^1,1. • • • : J] P r • • • . K,<) ft ~ ^ ~ E " " " > ^ 

2 Vi,i-Vi, 2 »=2 fc=l i=l 

V2.1-V2,* 

and, following the same steps as before, 

k k 2 

<n^>r +1 =E^ e n(w+E^-(^---^. 

j=l 2 Vi,i-Vi, a i=l i=l 

V 2 ,i-V 2 , z 

z 

P?(Vl,l, • ■ ■ V nA ) l[P°(V l!h ■ ■ ■ Vn,i) 
i=2 

k 

= ( TT Vrij )" + lower order terms 

3=1 
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